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. Abstract 

Q . The present paper is an extension of cond-mat/03 12673. 

CJ I The construction of a hybrid discrete-continuous model of layered su- 

perconductors is presented. The model bases on the classic Lawrence- 
Doniach scenario with admitting, however, long-range interactions be- 
tween atomic planes. Moreover, apart from Josephson couplings they 
' involve the proximity effects. The range of interactions can, in prin- 

■ ciple, be arbitrary large. The solutions corresponding to the range 

«) , K=2 are found. The mechanism of enhancement of superconductivity 

\C \ caused by the proximity effect and the presence of higher Joseph- 

' son couplings is shown. The physical meaning of coupling constants, 

. with particular attention paid to their sign, is discussed. For the case 

K=2 the interpretation in terms of microscopic interactions between 
Cooper pairs in different planes, as well as the relation to experimen- 
tally measurable quantities, such as the out-of plane effective mass 



Q , and bandwith, is given. 



1 Introduction 

Most of the high-temperature superconductors like e.g. YBCO or BiSCCO 
have a layered structure. Such a strongly anisotropic situation results in the 
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fact that the material properties and behaviour of the fields in the direction 
(say z axis) orthogonal to the layers are totally different from the behaviour 
in directions parallel to the layers. The simplest effective phenomenology of 
such systems is given by the 3D anisotropic Ginzburg-Landau model [2, 1, 3]. 

However, the abilities of the 3D continuum to describe discrete array of 
atomic layers are limited: for a layered superconductor such a description 
can be justified only when >> -s, where s is the (typical) interlayer spac- 
ing. The idea originally proposed by Lawrence and Doniach (LD) [6] is to 
consider the layered superconductor as an array of superconducting planes 
coupled by Josephson currents [6] flowing perpendicularly between adjacent 
layers. The planes themselves are described in terms of the 2D Ginzburg- 
Landau phenomenology. Under such circumstances the GL order parameter 
■0 appears as a function of two continuous variables (say x and y) and one 
discrete variable n - the index of the layer. The corresponding form of the 
free-energy functional was proposed in [6]. The model turned out very suc- 
cessful, and was subsequently modified and enriched in many papers and 
textbooks [7, 9, 8, 10]. A number of solutions for LD model are given in [12] 
or [11]. 

As compared with the 3D anisotropic GL, the LD model can be considered 

the other extreme case, justified when ^c<s. However, both the models 
together do not cover the entire range of the coherence lengths in relation to 
the spacing s: there is a gap between the models corresponding to coherence 
lengts which are too small for GL and too big for LD. 

The higher grade hybrid model, proposed in the present paper, patches 
the above gap by admitting direct couplings not only between adjacent, but 
also between more distant atomic planes. Such couplings can be particu- 
larly expected due to strong interplanar electronic correlations. The low- 
temperature correlation lengths of many superconducting materials, such as 
e.g. BiaSraCaCuaOg, with = 28A, or TlaBaaCaCuaOg, with = 23A, 
do not justify the restriction of interplanar couplings to nearest neighbours 
only. Since in both cases the interplanar spacing s is about 15A, the ra- 
tio (^c/s ranges, roughly speaking, from 1.5 to 2, and this suggests the need 
for interplanar couplings including nearest and next- nearest planes at least. 
The situation here is quite analogous to that in the theory of crystaUine lat- 
tices: in general the nearest neighbour couplings alone are not sufficient for 
adequate description of lattice dynamics [16]. 

The paper is organized as follows. In the next sections we give a formal 
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presentation of our hybrid model with couphngs among more distant neigh- 
bours - up to a (given, but arbirary) range with emphasis on = 2. 
Subsequently we determine the ground state solutions and discuss their sta- 
bility in function of admissible coupling parameters. For the sake of brevity, 
in other (frequently important) respects we keep our considerations as close 
to classic LD model as possible, avoiding e.g. hcterostructurcs, non-uniform 
intcrplanar spacings, or else gradient couphng [9]. Those factors will be dis- 
cussed elsewhere. 

In Section 7 we discuss the relation between 2D and 3D superconductivity 
in our model. The presence of interplanar couplings modifies the in-plane 
superconductivity. Depending on the coupling constants the modification can 
turn out supression or enhancement. In particular, even overcritical planes 
can, under appropriate couplings, result in subcritical superconductivity of 
the whole array. 

If the nearest-neighbour coupling dominates, one can pass to = 1. 
The resulting hybrid model with K = 1 offers, however, more possibilities 
than the classic LD; the latter obtains as a specialized case of this variant. 
The relevant questions, including the relation to other existing models, are 
addressed in the Section 8. The final Section 9 is devoted to discussion of 
the physical significance of the coupling constants, with particular attention 
paid to their signs. 

2 The hybrid model of grade K 

Let us regard the layered superconductor as a one-dimensional chain of 
atomic planes with Josephson's bonds between them. Such bonds will be 
called J-links. The interplanar distance will be denoted by s. We assume the 
following convention for indexing the planes and links. If we locate the point 
2; = at an atomic plane, then the 2;-coordinate of any plane, equal ns, may 
be represented by the integer n, while the z-coordinate of the center of any 
interplanar gap, equal Is, by the half-integer I. Choosing the point 2; = at 
the center of an interplanar gap - we index the planes by half-integers and 
the gaps by integers. 

Let us consider the free-energy functional T for a layered superconductor. 
We shall denote by ipn the order parameter associated to the layer indexed 
by the number n. Its complex conjugate (c.c.) will be denoted by The 
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symbol rriah and rric will denote the in-plane and out-of-plane effective mass of 
superconducting current carriers, respectively. We start from the free energy 
functional of the following form 

J^o + ^s + ^/B^d^x. (1) 

The term J^o describes the normal state, while the superconducting one. 
The supercoducting term is composed of two parts: 

J^s = :Fp + Tj, (2) 

where the part 

•^P = "5 H / dxdyFn (3) 

describes the contribution of atomic planes, while !Fj corresponds to intcrpla- 
nar Josephson's bonds. For any plane indexed by n the free energy density 
Fn has the 2D Ginzburg-Landau form (in general, the parameters can depend 
on n) 

F„ = + ao\^n\' + Ipi^nW (4) 

2mab 2 

where we have introduced the 2-dimensional continuous operator D (covari- 

ant derivative) 

ic* 

Dp^dp- -^Ap, P^x,y. (5) 

The form (4) of the functional J^p already ensures its invariance with 
respect to the gauge transformation 

A ^ A' = A + VA, 

V'n ^ < = V'ne^^^, (6) 

(2-dimensional A and V for this case). The standard variational treatment 
of the functional JF^ with respect to Ap, p — x,y, gives the standard in-plane 
components of the supercoducting current 

ic*ti — — 6*^ 

jp = -7, {i^ndp'^n - tpndptpn) ^plV'nP, P^X,y. (7) 

2mab mabC 



Let us now construct the term jFj in the free-energy. In general it is 
a functional which can depend on all ipn,4'n and on the vector potential A. 
The simplest gauge invariant expression for the energy of J-link between n-th 
and {n + q)-th planes will have the following form 

^qn = ^{Cqni'ni'n + Tjgni'n+qi'n+q " {lqn1pn1pn+qe~'^''" + CC.)}, (8) 

where the exponent Pgn is defined by the formula 

g* r-(n+q)s 

Pqn= T- ^zdz. (9) 

liC J ns 

Let us note that for the model invariant with respect to the time reversal we 

have 7q„ = 7g„. In general, the coupling parameters (, 7], 7 as well as in-plane 
parameters ao and f3 can be different for different planes (which is the case 
for superconductors composed of various atomic planes). However, in this 
paper we shall consider the array of identical planes, hence the parameters 
will not depend on the index n. That implies rjq — Q. Hence instead of (8) 
we shall use 

^qn = ^{Cgd^nP + l^n+gH " 7q (^nV'n+.e^^P^" + CC.)}. (10) 

Let P denote the (finite or infinite, but ordered) set of all indices of planes, 
and let Q = {1, 2, K}. The planes connected with the n-plane by Joseph- 
son coupling select the following subset of Q 

Pn = {qeQ: (n + g)eP}. (11) 

Thus, the gauge invariant functional J-j for the hybrid model of grade K has 
the following form 

y^j^s^Yl / dxdysgn, (12) 

neP qePn ^ 

with Sqn given by (10). The coupling parameters C,q and 7^ vanish ior q> K. 
Every J-link is represented in (12) by exactly one term. 



3 Comparison with anisotropic GL model 

To compare our hybrid model (HM) with the continuum GL model, we shall 
consider the infinite medium. In that case the summation index may be 
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shifted by the integer q. Moreover Pn = Q. This imphes that 

Fj = ^ EE{2(C, - 7qMn\' + 7,|^n+.e-'^- - (13) 

n q 

For very weak field and very slow dependence of ipn on n we have the 
correspondence rules which, in the long wave limit, allow us to pass from the 
hybrid to the 3D continuum: 

' En ^ i / t^^, 

^(V'n+9-V'n)^V';(^), (14) 

Applying the rules to the functional (2) with (4) and (13) one obtains 

:fs^I ^M^iDV'r+«oiV'r+^/?iV'r+E[(c.-7,)iV'r+^g^.^,p3V'r]}- 

(15) 

Hence, we have the following relation between rUc - the effective mass in 
2;— direction (in the anisotropic GL model) and the coupling parameters 7^: 



1 

rric h a 



(16) 



The presence of J-links modifies also the parameter ao to the form: 

a = ao + E(Cg-7q)- (17) 



The condition of weak dependence of ipn on n, valid under usual assump- 
tion of nearly uniform states, is crucial for the transition from discrete to 
continuum picture in the z-direction. For strongly oscillating dependence, 
in particular for alternating ground-state solutions discussed in the sequel, 
transition (14) can no longer be estabhshed. Nevertheless, sufficiently weak 
and sufiiciently smooth excitations from a given ground state can again be 
described by an effective 3D anisotropic GL model, specific for this ground 
state. The equations (16) and (17) remain valid only for excitations with 
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respect to the uniform ground-state. In the case of phase-oscillatory ground 

o 

states ipj^, making use of the substitution 

(18) 

with a suitable rcintcrpretation of the correspondence rules (14), one arrives 
at the modified equations, which can be obtained from (16) and (17) through 
the replacement 

7q ^ 7qC0s(ge). (19) 

The requirement of the overall stability of the system is responsible for the 
positive sign of the effective mass rric- 



4 Field equations 

By computing the variation of the functional JF with respect to one obtains 
the Maxwell equation for the 2;— components of current density and curli? 

-J{z) = ^{cmlH),, (20) 
c 47r 

where 

se* - 

J(^) = ^ E E Mn^n+qe-'P"- - C.C.}x,n(^), (21) 

neP qePr, 

Pn is given by (11), and the quantity Pgn by (9). The symbol Xqn{z) denotes 
the characteristic function of the interval [ns, {n + q)s\. Let us note that 
for any layer ns < z < {n -\- l)s the exppression J{z) does not depend on 
the value z; only the ends of the interval are important. To better see the 
structure, let us first extend the set Q on the negative values 

g = {-K,...,-2,-l,l,2,...,i^}, (22) 

and introduce the symbols 

_ f 1 if (n + g) eP , , . 

1 otherwise. ^ > 
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Then the expression for Joscphson current J; describing tunnehng across the 
interplanar gap indexed by / (half integer if P contains integers) will have 
the form 

sc* - 

Jl=W^Y.Y. ilgf^qn^n^n+qe-'P"- - C.C.}, (24) 
^^'^ geQnePi^ 

where 

Piq^ {neP -.TKl <n + q}. (25) 

By computing the variation of the functional J-g with respect to one 
obtains the equations 

1 

-7; ^^i^n + + P\^n\^^n- 7;12%'^m^n+qe~'P'^" = (26) 

2mab 2 
where instead of ckq we have introduced 

a ^ ao+^Yl (^qnCq (27) 
which, for finite P, depends of n. 



5 The ground states 

Let us now consider the plane-uniform states of HM in the absence of mag- 
netic field. The order parameter is then independent of the in-plane variables 
and the net supercurrents vanish. In detailed calculations we shall focus our 
attention on the grade K — 2, which seems to be sufficiently illustrative to 
grasp the idea on what is going on. The generalization to grades K > 2, 
although more complicated algebraically, is straightforward. In the region 
far from the boundary all the coefficients aqn = 1. For K = 2 the condition 
of vanishing Josephson current is equivalent to 

7l('0n'0n+l - C.C.) + 72('0n'0n+2 + V'n-lV'n+l " C.C.) = 0, (28) 

and the equations (26) take the form 

aipn + /3\lpn\'^1pn " ^[7l(V'n+l + V'n-l) + 72(V'n+2 + lpn-2)] = 0. (29) 
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Wc shall look for solutions with constant amplitude and difference of phase 
between adjacent atomic planes, so we use the ansatz 

V'n = Ce^"^ (30) 

The result is the equation 

a + /?C2 - 7i cos e - 72 cos 29 = 0, (31) 

with the condition 

7isin^ + 272sin2^ = 0. (32) 

Solving (32) with respect to 9, we obtain 3 variants: (al) ^ = 0, (a2) 9 — 
and (a3) cos^ = — 7i/472. The solution C to (31) has the form = —a*/P, 
where a* depends on the variant. The variant (al) implies the uniform 
solution to (29): 

= C, (33) 

with 

a* = "0 + Ci + C2 - O'l - 72, (34) 

what is the case of equation (17) for grade K=2. In the variant (a2) the 
solution to (29) has the alternating form 

= ±C, (35) 

with 

a* = ao + Ci + C2 + 71 - 72- (36) 

Finally, in the variant (a3), the solution exists if the parameters 71 and 72 
fulfill the relation 

|7i|<4|72|. (37) 

Then the parameter a* is connected with the coupling constants by the 
formula: ^ 

a* = ao + Ci + C2 + 72(1 + A)- (38) 

072 

There are two independent solutions 

^jJ^^Ce^'^, e = arccos(-^). (39) 

472 
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They will be referred to as the phase modulated states. The solutions de- 
generate at the extremities |7i| = 4| 72I. 

So far we confined our discussion to the existence of solutions which could 
serve as candidates for the ground state. However, unstable candidates, even 
uniform ones, have to be rejected. Note that alternating states were already 
found in the literature on different grounds - phcnomcnological [7] as well 
as microscopic ones [8]. The question of stability of the solutions will be 
addressed in the next section. 

Let us note that the condition sin^ = admits the solutions (33) and 
(35) for hybrid model of any grade K. Contrary to that, the analogues of 
the conditions (a3) and (37) can deliver, depending on the grade and on the 
coupling parameters, any number from to 2[K — 1) modulated solutions. 



To examine the stability of the solutions presented in Section 5, we shall 
analyze the Hessian matrix of the free energy J^s describing small deviations 
from the ground state. The problem reduces to examining the function 



The solutions found in the previous section arc stationary points of 8. If 
72 = 0, then the stability of the solutions depends on the sign of 71. If 
72 7^ then the sign of the respective eigenvalue depends on the values of 
the parameters 71 and 72 . For 9 different from and tt the dependence has 
form according to the function 



The straight lines 71 + 472 = and 71 — 472 = divide the plane 71, 72 into 
4 regions (see Fig. 16). 

As explained above, both the uniform and the alternating solutions always 
exist. However, in the region 71 > 0, 71 + 472 > 0, only the uniform solution 
(33) is stable, while in the region 71 < 0, 71 — 472 < 0, only the alternating 
solution (35) is stable. The region (S): 72 < 0, 472 < 71 < —472, excludes the 
stability of both the uniform and the alternating solutions but, in contrast 
to that, ensures the existence and stability of the modulated solutions (39). 



6 Stability 




(40) 



f {11,12) = 72(71 + 472) (71 - 472). 



(41) 
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In the region (N): 72 < 0, —472 < 71 < 472, the modulated solutions exist 
but are unstable. 

Let us note that in the regions (NW) and (W) of stabihty of the alter- 
nating solution (35) one can apply the construction of bi-continuum solution 
presented in [13]. Note also that, irrespectively of the values and signs of 
7i, 72, a stable ground state solution always exists. 

7 Enhancement of the superconductivity 

The asociation of the formulae (34), (36) and (38) with the regions of stabihty 
of the ground state shows that, for suitable relations between the couphng 
constants 71 and 72, one can make the parameter a* more negative than 
ttQ. In consequence, the 3D superconductivity can turn out enhanced with 
respect to the 2D superconductivity in the atomic planes. Such a possibility 
has been indicated in the literature [4]. The presence of coupling constants 
C,q allows to describe the proximity effect between atomic planes. 

The corresponding shift of the parameter a depends on both C's and 7's. 
While the contribution from (^'s reduces always to simple renormalization 
of tto by the same additive term C,i + ^2, the contribution from 7's in (34), 
(36) and (38) varies. To expose the dependence on 7's, we shall count the 
enhancement with respect to the origin located at a given by (27). It is 
convenient to introduce the polar coordinates 7 and </? in the plane 71, 72: 

7i = 7COS99, 72 = 7sin(/?, (42) 

and the notation ^ 

</7o = arctan(-). (43) 

The quantity Aa — a* — a as a, function of the coupling angle </? is plotted 
in Fig. 16 (the numerical values are computed for 7 = 1). 
In the uniform state (34) we have 

3 

Aa = \/27 sin(93 — -tt), —(fio^'P^T^ — ^o- (44) 

The minimum Aa (hence the maximum enhancement) is reached at </? = 7r/4 
and equals —y/2j. 
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In the alternating state (36), in turn, one obtains 

3 

Aa = -\/27sin(v? + -tt), v'o < tt + <^o (45) 

with the minimum value —\/2^ reached at 99 = |7r. 

The enhancement for the phase modulated state (38) is represented by 

Act = 7(1 + ictg^^j) sin^?, -'K + ipQ<ip< - LpQ. (46) 

o 

In this case the minimum equals —7 and is reached at 99 = — 7r/2. Hence, 
the maximum enhancement in this case is smaller than in the uniform and 

alternating states. 

Let us note that, in the enhancement mechanism discussed above, the 
2D superconductivity of the planes is not a prerequsite for the 3D supercon- 
ductivity of the array of the planes. In fact, one can obtain the negative 
a* starting from positive a^. This is in concordance with ideas expressed in 
Anderson's discussion of his Dogma V in [5]. 



8 Special cases 

The HM with K=l has in general two coupling parameters C, and 7. If they 
are equal to one another, one obtains either the Lawrence-Doniach [6] or 
Theodorakis [7] model, depending on the sign of 7. The formula for the 
effective mass in z-direction simplifies to 

--U^V (47) 

Due to 7 > the LD model has only one stable solution, namely the uniform 
ground state. The effective parameter a = ao, and the interplanar coupling 
gives neither enhancement nor suppression of the critical temperature. Al- 
though the Josephson current coupling places the model on the enhancement 
side, the efect is precisely annuled by C = 7- 

Where the negative 7 is concerned, the ground state solution is also 
unique. Contrary, however, to the LD case, this solution is alternating, the 
proximity effect is present, and the effective parameter a — ao — 2\^\ results 
in enhancement of the supeconductivity. 
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9 Discussion 



In the present section we shall discuss the physical meaning of the couphng 

constants 71, 72 with special attention paid to their signs. 

1. First of all it is expedient to note a certain symmetry between hybrid 
systems with negative and positive 71. To that end consider the transforma- 
tion 

r = (-l)-V'n, 

< A' = A, (48) 

. = (-1)%, 

with all the remaining quantities kept fixed. This transformation switches 
between systems with all the constants ■jq transformed according to the parity 
of q: the constants remain identical for even g's and change sign for odd g's; in 
particular 7^ = —71 and 72 = 72- At the same time the order parameter ip^ is 
transformed by the sign- alternating factor; in particular the states of uniform 
order parameter are transformed into alternating ones, and vice versa. 

An inspection of the equation (8) shows that the energy of the Josephson's 
coupling is invariant with respect to transformation (48). This invariance 
extends to the whole free energy functional (1). Note, in particular, the 
correspondig symmetry in Fig. 1 and Fig. 2. Two hybrid systems related 
by the transformation (48) have identical ground state energies, excitation 
energy spectra, and all the thermodynamical properties, including Tc. 

2. The second coupling constant 72 represents the microscopic interac- 
tions between Cooper pairs located in next-nearest layers. Similarly as in the 
case of 7i, corresponding to ^3 in [8], the sign of 72 reflects the attractiveness 
or repulsiveness of those interactions, both signs being equally well admissi- 
ble on physical grounds. Contrary, however, to 71, the sign of 72 can not be 
relativized by the transformation (48). To discuss this situation let us first 
select one of the admissible ground states determined in Section 5, and then 
consider the related excited states resulting from ^-dependent variations of 
the order parameter with respect to the selected groumd state. 

3. For simplicity we start here from the case of sufficiently weak coupling 
between next-nearest layers, so that I72I < |7i|/4, what corresponds to the 
western and eastern sectors in Fig. 1, with the alternating and uniform 
ground states, respectively. The strong coupling sectors will be commented 
in the sequel. 
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Understood as the excess energy over the ground state and considered as 
a fuction of the wave vector, the excitation energy equals 



efc = (A;^ + /c^) + |7i|(l-cosM +72(1- cos 2A;^s). (49) 

Note that the contribution from 71 is always positive. 

In the present work it will be sufficient to consider only excitations with 
— ky = 0. They are characterized by the effective mass 

mc=-(|7i|+472)-' (50) 
s 

in the vicinity of the ground state A; = 0, and by the bandwidth 

W^ = 2|7i|; (51) 

the latter furnishes an alternative interpretation of 71 in this case. Provided 
that the bandwidth W and the effective mass rric (as well as the interlayer 
distance s) are known from measurements, the coupling constant 72 can be 
determined from the equation 

T. = i(^-H') (62) 

where 

2s^mc 

M = (53) 

represents the effective mass, scaled conveniently for the sequel. The inequal- 
ities 

2 < MW < 00 (54) 

hold in the western and eastern sectors. 

The quantity Wld = 4:/M represents the bandwidth for 72 = 0, which 
corresponds to the LD model. This bandwidth is entirely determined by 
the spacing s and the effective mass rric alone, what is one of artifacts due 
to oversimplification of the model resulting from 72 = 0. In our model the 
effective mass and the bandwidth are independent measurable parameters. 
According to equation (52), the sign of 72 is directly determined by the rela- 
tion between the measured bandwidth W and the parameter Wld calculated 
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from the effective mass: when W < Wld, the constant 72 is positive, and 
when W > Wld, it is negative. 

4. Now we shall briefly discuss the case of strong next-nearest coupling, 
defined by the relation (37), represented by the remaining two sectors in Fig. 
1. The dispersion (49) is no longer monotonic, and in place of the simple 
expression (51) for the bandwidth we have 

|7i| + 2|72| + -^, (55) 

5|72| 

the formula common for both sectors. However, the ground-state effective 
masses in the northern (72 > 0) and the southern (72 < 0) sector differ from 
one another. 

In the northern sector we have, according to the results of Section 5, two 
stable solutions with generically different energies: the ground state and the 
metastable state. For the ground state effective mass the equation (50) is 
still vahd. Combined together, the equations (55) and (50) allow to estabhsh 
the inequalities 

1 < MW < 2 (56) 

which restrict the admissible range of {rric, W) in the northern sector. Pro- 
vided that these restrictions are satisfied, the equations (55) and (50) furnish 
the unique solution for the coupling constants 

The southern sector is characterized by two symmetric phase- modulated 
ground states. The corresponding effective mass is given by 

s2 I67I - 7f 

and in place of the inequalities (56) we obtain 

1 < MW < 00. (60) 
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Finally, under the above restriction, the coupling constants are uniquely de- 
termined from (M, W) as 



5. For a known sector, making use of the relevant pair of equations (51)- 
(52) or (57)-(58) or else (61)-(62), one can uniquely determine |7i| and 72 
from experimental values of the effective mass M and the bandwidth W. 
However, identification of the sector with the aid of inequalities (54), (56), 
and (60) is possible only partially: eg. west, east and northern sectors can 
be mixed with the southern one, while MW = 2 constitutes the border 
between west and east sectors on one side, and the northern sector on the 
other side. Nevertheless, more detailed knowledge concerning identification 
of the sectors can be gained from closer examination of the excitation spectra. 
In particular, a feature which distiguishes sectors N+S from W+E is the 
presence of a strong in-band singularity of the density of states N{ek). 

6. In consequence of the symmetry (48), one can not experimentally dis- 
criminate between sectors W and E by measurements based on the ground 
state or excitation energies, or else on any physical quantities determined 
by the energy spectra. Nevertheless, such a discrimination should be possi- 
ble with the aid of suitably chosen measurable physical quantities which are 
not invariant under (48). A similar problem has already been addressed by 
Lazarides, Schneider, and S0rensen in [8], where solutions differing by the 
sign-alternatig factor and yet having identical energies have been found in 
a BCS-like microscopic model of multilayer superconductors, depending on 
the sign of the coupling between adjacent layers. Attractive interaction leads 
to positive coupling constant and uniform ground state, while repulsive in- 
teraction results in negative coupling constant and alternating ground state. 
The authors suggest two experimental settings which can determine the neg- 
ative sign of this constant: a superconducting loop able to detect trapped 
half-integer flux quanta, and a DC SQUID conflguration able to detect the 
presence of 7r-junctions, corresponding to alternating solutions. 




(61) 




(62) 
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Figure 1. The regions of stability: uniform 0, 
alternating @, and phase modulated @ states. 
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Figure 2. Enhancement of superconductivity 
by Josephson couplings: Aa vs. </?. 
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